Evolution of Kerr-Schild type initial data for binary black holes using the horizon 

penetrating Teukolsky equation 




I. INTRODUCTION 



There is a lot of current interest in studying the colhsion of two black holes since these events could be primary 
sources of gravitational waves for interferometric gravitational wave detectors currently under construction. Numerical 
simulations are considered as the key to model such collisions. Unfortunately, current computer limitations cause these 
simulations to be short lived, so to study the actual merger of two holes, one needs to have initial data that represents 
two black holes close to each other. 

Early attempts to provide binary black hole initial data solved the initial value problem of General Relativity using 
conformally flat slices [Q. This restriction was too strong; it was unable to incorporate spinning holes. Lately, some 
effort has been geared towards obtaining initial data based on the Kerr-Schild form of black hole solutions. The main 
advantage of this form is that the spatial slices of these solutions are horizon penetrating, which make them ideal 
for the application of the "excision" method for evolving black holes. This form of the solution also naturally 
incorporates boosted and spinning black holes. 

There is also another way Q to obtain initial data for a black hole binary, without solving the initial value equations 
directly. One could, in principle, consider a scenario with matter fields coupled to gravity such that there are no black 
holes to begin with. One could distribute the matter fields in such a way that their gravitational effects are like those 
of two stars that are separated far apart but each is critically close to gravitational collapse. It is quite possible that 
such a simulation would not be plagued with the subtle problems that ocurr when one attempts to evolve two black 
holes. One could set the parameters of the simulation in such a way that when the two stars come close to each other, 
they collapse upon themselves and form two black holes. This final slice could be used as initial data for a binary 
black hole code. 

However in this paper, we consider the Bishop et al. j2j family of initial data wherein two non spinning Kerr-Schild 
black holes are superposed. This proposal has the property that in the "close limit" in which the separation of the 
holes is small, the metric is given by a distorted Kerr-Schild black hole. The evolution will be carried out using a 
recently introduced horizon penetrating Teukolsky formalism [Q . We present Teukolsky waveforms from this evolution 
and the amount of energy radiated away by the system. The results we obtain are in agreement with those obtained 
by Sarbach et al. ^ from their evolution of the same initial data family using a Kerr-Schild type Zerilli equation. 



For the sake of completion, in the first subsection of this section we reproduce the calculations followed by Sarbach 
et al. ^ for an exact solution of the linearized initial value problem. 
The KS space-time metric is defined by 



where is a null vector. The rj^^, is taken to be the Minkowski metric with coordinates {t,x) = (i, x*) such that 
r]tt — —1, riti = and rjij — 6ij. The null vector fc^ satisfies kt — —1 and k''ki = 1, where fc' = S^^kj. The three metric 
and extrinsic curvature with respect to a slice t = const, are 



II. KERR-SCHILD INITIAL DATA 



- 2Vk^k, 



(1) 




(2) 



where a — {\ — 2V) is the lapse and where V refers to the flat metric Si 
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Bishop et al.'s solution procedure |^ consists in inserting equations |] and ^ into the constraint equations and to 
solve the resulting equations for V, V — dtV and ki — dtki, where ki is assumed to be given. 

A. Two black hole data 

Here we review the discussion of equations |l| and || as followed by Sarbach et al. where a particular ansatz is made 
for ki representing two nearby non-rotating and non-spinning black holes. A single Schwarzschild black hole can be 
written as 

with = 1/r. For two black holes, Bishop et al. choose 

n^) = 1 r + 



where and X2 are the positions of the black holes that have mass Mi and M2. If the black holes are located at 
Xi = (0, 0, ai) and X2 — (0, 0, a2), with ai > > 02, may be expanded in a sum over multipoles: 

0(r, I?) = 2^ ^ P£(cos d), (3) 

where Pi denote standard Legendre polynomials and where (r, 1!), ip) are polar coordinates for x. Note that the 
expansion (j^) is only valid for r > max{ai, — a2}. We define the separation parameter as 

ai — 02 



M ' 

where M = Mi + M2 is the total mass. Now, using the center of mass condition MiOi -I- M2a2 ~ 0, the close limit of 
becomes 

(/. = — -f -^P2(cosi9) -f 0(eV»- )■ 

As a result, to first order in e^, ki is given by 

fc, = -l, fc^=e^^v^p,, (4) 
r 

where here and in the following, A = "iJ, (ys. The remaining amplitudes are expanded according to 

M ■ ... 

+e'v{r)P2, V^e^v{r)P2, k a ^ e'' k{r)V aP2 ■ (5) 

r 

Plugging this into the constraint equations, and keeping only linear terms (of the order e^), one obtains the equations 



3M / 2M\ • 3 GMM1M2 
-v + —[l + U - - 



, 4 ^&M\ 6MM1M2 , 
Q = -v v^ —k (r-M), (7) 

2M; 6MM1M2 

= -Mk' + V + —k + (8) 

Here, a prime denotes differentiation with respect to r. The system can be re-expressed as a single second order 
equation. Introducing the dimensionless quantities x = r/M and ^ = M1M2 / M'^ , this equation reads 

Q = -v^^--v^ + ^v + ^{ix + 2). (9) 
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Once we have solved this equation, the remaining amphtudes k and v are obtained from (||) and respectively. 
A particular solution of is given by 




In order to find the solutions of the homogeneous equation, one performs the transformations x — v{x) — z^u{z), 

which yields the Bessel differential equation 

~ Z^Uzz + ZUz — (16 + z^)u. 

The solutions arc a linear combination of the Bessel functions 14(2:) and Ki{z). Thus, the general solution to (^) is 
v(x) — fiv{x), with 

Hx) = -|fl-- + 4 + 4 + ^K.iVM/x)) + ^/4(V247^). (10) 

o \ X X X X J X 

Choices for A\ and Ai are discussed in the following section. 



B. Initial data for the Teukolsky function 



The initial data for the metric and extrinsic curvature obtained above will be evolved using the KS type, penetrating 
Teukolsky equation. This equation was derived by Campanelli et al. in order to study gravitational perturbations 
close to the horizon. The equation, along with the associated geometrical quantity, the Newman-Penrose complex 
null-tetrad is regular at the horizon, and therefore has the capability of evolving through it. Such an equation is ideal 
to study the effects of "excision" as used by full numerical, binary black hole collision codes, like MAYA [||. It should 
be noted, that the penetrating Teukolsky equation is based on a rescaling of the transformed Kinnersley tertrad (to 
KS coordinates). This rescaling enables the tetrad to be regular at the horizon. Therefore, the penetrating Teukolsky 
equation is not just the original Teukolsky equation in KS coordinates. 

To obtain initial data for the Teukolsky function, we use the formula developed by Campanelli et al. ||] in terms of 
the perturbative three metric and extrinsic curvature as shown below. "04 is a Newman-Penrose scalar, that captures 
information about outgoing gravitational radiation, and it is related to the Teukolsky function as, ip = r'^'ip4 (for 
a non spinning black hole). 



-04 



K 



4 [Rji - KjpKf + KKji]^^^ n[°m^ln[^ 

doRji - 2K^idoKj)p - 2N(Q)j\jp 



(1) 



-|-8n[°TO-'ln''m' 



(1) 



(11) 



(12) 



where 



i'jkl 



2K,[kKi],] Tn^m'^n' + 4 



K 



(1) 



(13) 



'^We recently learnt that Carsten Kollein at Albert Einstein Institute performed a perturbative evolution of this same family of 
initial data in his M.Sc. thesis (unpublished), using the Teukolsky formalism. However, the Teukolsky equation used by Carsten 
was the Boyer-Lindquist Teukolsky equation transformed to KS coordinates, and not the penetrating Teukolsky equation. We 
thank Jorge Pullin and Manuela Campanelli for providing us with a copy of Carsten's work. 
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FIG. 1. The initial Teukolsky function, as function of r and 6. The value of the constant A2 here is 0. 



^(0) = (fi'fc'err) "'^^^ the zcrotli order lapse, in these equations should be taken to be related to that of the 
original penetrating tetrad as — ^(0) "-orig ' = "-orig + '•^■'n^. Latin indices run from 1 to 3, and the brackets 
are computed to only first order (zeroth order excluded). The derivatives involved in the above expressions can be 
computed in terms of the initial data on the Cauchy hypersurface as, 

aoif = iV(o)ifp,ifP«-V2iV(o), (14) 
doR = 2KPidoKpq + AN(o)KpqKPK"i - 2KdoK, (15) 

doRijki = -47V(o) ^K,[kRi]j - Kj[kRi]x - ^R {K,[k9i]j - Kj[k9i\i) | (16) 
+'^9i[kdoRi]j - 2gj[f,doRi]i ~ 9i[k9i]jdoR + '^Ki^^doK^j - 2Kj[j,doK^^, 

and, 



doKij = iV(o) 



% + KK,, - 2K,pKPj ~ 7V(-iV,V,7V(o) 



(1) 



(17) 



We calculate the Teukolsky function and its time derivative ■0, on the initial time slice, using the above formulae 
on expressions for the metric and extrinsic curvature that we produced before. The algebraic expressions for the 
initial data of the Teukolsky function are too large to include in this paper. Instead, we include some surface plots. 

We still need to discuss the choices for the unknown constants Ai and A2 that appear in equation |l^ and therefore 
in expressions for tp and ^p. It turns out that in order for the Teukolsky function to have a proper fall-off behavior 
for large r, Ai must equal —12. This choice for Ai also ensures that v(r) vanishes at spatial infinity. A2 on the other 
hand can be chosen at will. In their paper, Sarbach et al. provide a very clear explanation regarding the meaning of 
A2 based on the original assertion by Bishop et al. They show that the meaning of A2 is related to the position of 
the apparent horizon, through a genuine choice (not just a gauge choice) in the initial data. 

In what follows, we shall choose Ai = —12 and use multiple values of A2. In particular, we shall use A2 — and 
A2 = 10 for our results. 



III. EVOLUTION 



Like we mentioned in the previous section, the initial data for the Teukolsky function obtained above, shall be 
evolved using the penetrating Teukolsky equation. This equation is best suited for comparisons with current KS type 
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FIG. 2. The initial time derivative of the Teukolsky function, as function of r and 9. The value of the constant A2 here is 0. 

full numerical binary black hole codes like MAYA, because its coefficients are regular at the horizon. This equation 
was derived by Campanelli et al. [Q, and we reproduce it here (with the Kerr parameter, a — Q), 



To implement the penetrating Teukolsky equation numerically, we reduced it down into a 2+1 dimensional one, using 
a decomposition of the Teukolsky function into angular modes, = S^/'me™'^. We used Lax-Wendroff technique to 
numerically implement these simplified set of equations. 

Results of evolutions for two equal mass KS holes are shown in the included figures. The inner boundary for the 
evolution was chosen to be inside the horizon, at 1.5M. In figures || and ^ we show the Teukolsky function as a function 
of time, extracted at r = 50M. We show waveforms for two values of A2. Quasi-normal ringing is self-evident in these 
figures as expected. We also calculated the energy radiated for these evolutions, and we get amounts in agreement 
with those obtained by Sarbach et al. 



If we pick A2 such that the energy radiated is at a mininum (this occurs at A2 = 7.87), and we pick the conformally 
flat separation between the holes to be about 0.9M, the value we obtain is about 10~^M, which is consistent with 
the value obtained from the evolution of Misner data (with the same separation L — 0.9M or /xq = 1.5) using the 
original Teukolsky or Zerilli frameworks [pj. 



We evolved the initial data family of Bishop et al. in the close-limit by treating the spacetime as a single distorted 
Kerr-Schild black hole. We performed the evolution using the penetrating Teukolsky equation. 

These results shall be used in the calibration of numerical codes based on the Kerr-Schild coordinates, e.g. MAYA 
1^. They are also in agreement with those obtained by Sarbach et al. from their evolution of the same initial data 
family using a KS type Zerilli equation. 




(18) 



(4r + 6M) + 2(2 cot^ d - 1)V^ 



0. 



E = e^M (1.31 X lO"'^ + 1.22 x IQ-^Al - 1.92 x IQ-^Az) 



IV. CONCLUSIONS 
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FIG. 3. The radiated waveforms at r = 50M for A2 = 0. As usual in close limit collisions, the waveform is dominated by 
the fundamental quasi-normal mode. 
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FIG. 4. The radiated waveforms at r = 50M for A2 = 10. As usual in close limit collisions, the waveform is dominated by 
the fundamental quasi-normal mode. 
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